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Molecular and atomic gas in the young TeV 
7-ray SNRs RX J1713.7-3946 and RX 
J0852.0— 4622; evidence for the hadronic 
production of 7-rays 



Yasuo Fukui 



Abstract The interstellar molecular clouds are the site of star formation and also 
.^ i the target for the cosmic ray protons to produce 7-rays via the hadronic process. 

The interstellar atomic gas is enveloping the molecular clouds and may also be 
dense enough to affect the 7-ray production. In this Chapter, some of the basic prop- 
,S^ ' erties of the interstellar gas both in molecular and atomic forms will be reviewed. 

O-i. Then, it is presented that two young TeV 7-ray SNRs, RX J1713.7-3946 and RX 

Q I J0852. 0—4622, show good spatial correspondence between the 7-rays and the inter- 

H ■ stellar protons. The good spatial correspondence provides a support for the hadronic 

c/2 ' origin of the 7-rays in these SNRs. It is emphasized that both molecular and atomic 

hydrogen plays a role as targets for cosmic ray (CR) protons. The clumpy distribu- 
tion of the target interstellar medium (ISM) protons is crucial in the interaction of 
the supernova shocks with the ISM, whereas models with uniform ISM distribution 
are not viable. Finally, it is suggested that the dense atomic gas without molecules 
vQ . may occupy the dominant part of the dark gas in the local ISM. 

(N 



1 Introduction 



Recent progress in 7-ray imaging of supernova remnants (SNRs) is remarkable. In 
particular, several 7-ray telescopes, H.E.S.S., VERITAS, MAGIC, Fermi and AG- 
ILE, are obtaining high quality 7-ray images of high-energy objects including SNRs 
and PWNs in the Galaxy in the last decade (e.g., Gast et al. 2012; Hui 2010; Zanin 
C^ . 2009; Castro & Slane 2010; Giuliani et al. 2011). It was originally suggested in 

the middle of the 20th century that we are able to observe the 7-rays produced by 
the hadronic interaction between the cosmic ray (CR) protons and the interstellar 
protons which produce neutral pions decaying into 7-rays (Hayakawa 1956). 
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Fig. 1 NANTEN '-CO(7=1-0) integrated intensity distribution toward the W28 region (linear 
scale in K km s^ ' ) with TeV y-ray significance contours overlaid (green) 4, 5, 6c7 levels (Aharonian 
et al. 2008 Figure 2). The velocity ranges are 0-10 km s^' {!eft panel) and 10-20 km s^' {right 
panel). The radio boundary of W28, GRO J1801 -2320 and the location of the Hll region W28A2 
(white stars) are shown in the in the panels. 



The history of observations of the interstellar medium (ISM) spans over more 
than a half century and observations of the ISM have been made mainly at radio, 
infrared and optical wavelengths. The Hi observations at 21 cm and the CO obser- 
vations at 2.6 mm provide powerful tools to probe the major neutral component of 
the ISM. In particular, early CO surveys allowed us to make a comparison between 
the 7-rays and the molecular protons (Dame et al. 2001). The COS-B satellite ob- 
tained the 7-ray distribution of the whole sky in the 100 MeV-10 GeV range (Mayer 
et al. 1980). The correlation between the 7-rays and CO is good at a degree-scale 
resolution and was used to derive a X factor (=A^(H2)/W('^CO)) that relates the CO 
intensity to the proton density (Lebrun et al. 1983; Bloemen et al. 1986; Strong et 
al. 1988). 

The TeV 7-ray astronomy with ground-based atmospheric Cerenkov telescopes 
including H.E.S.S. opened a new era by detecting the highest energy 7-rays in the 
Galaxy. The TeV 7-rays may be produced by the CR protons whose energy is close 
to the knee and have a potential to probe the highest energy Galactic CRs and the ac- 
celeration site. The H.E.S.S. survey of the Galactic plane which detected more than 
50 7-ray sources was compared extensively with the NANTEN CO observations and 
the region of the W28 SNR was found to show the best spatial correspondence with 
the TeV 7-rays as shown in Figure 1 (Aharonian et al. 2008). It however remained 
as a puzzle why the other 7-ray sources do not show clear correspondence with the 
CO. For instance, the most remarkable TeV 7-ray SNRs, RX J1713. 7-3946 and 
RX J0852. 0—4622, were not identified conclusively to be associated with CO (e.g., 
Aharonian et al. 2007b). In this Chapter, we show that this puzzle has been solved 
by a careful analysis of the ISM protons including atomic protons. 
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2 Interstellar medium 

The interstellar medium consists of gas and dust. The mass ratio of the gas and dust 
is 100:1. The gas is either neutral or ionized. The ionized gas is a minor component 
in mass, since recombination is rapid in the Galactic disk for the general ionizing 
UV radiation field from OB stars. Therefore, the neutral gas is the dominant com- 
ponent in mass. The neutral gas is mainly either atomic or molecular hydrogen. 
Because of the higher cooling rate at higher densities, the denser gas tends to be 
cooler; the cooling rate is proportional to density squared, whereas the heating rate 
linearly depends on density. 

The atomic gas (Hi gas) has a density from 0.1 cm^^ to 100 cm^^ with kinetic 
temperature from 30 K to 3000 K. It is generally believed that such interstellar 
medium consists of two dynamically stable phases, either "high temperature and low 
density" and "low temperature and high density" (Spitzer 1978). It is yet possible 
that a significant amount of transient Hi gas at an intermediate density range exits, 
since the ISM gas is highly turbulent and far from the dynamically equilibrium (e.g., 
Koyama & Inutsuka 2002). It is generally assumed that the Hi emission is optically 
thin as suggested by the non-flat Hi line profiles (Dickey & Lockman 1990), whereas 
it is possible that there exits opaque Hi emission from dense Hi. Such Hi gas may 
be important in a shell compressed by the stellar winds in SNRs, although the nature 
of such Hi gas has not been discussed often in the literature. 

The molecular gas has density greater than 10^ cm^^ with typical temperature 
of 10 K (e.g.. Goldsmith 1987). The low temperature is in part due to the efficient 
molecular cooling and the strong shielding of stellar photons by dust grains. The 
major heating of the molecular clouds is made by CR protons of a 100 MeV range. 
The main constituent of the molecular clouds is H2 but H2 does not have a radia- 
tive transition observable at low temperature of 10 K. We use usually the rotational 
transition of CO at 2.6 mm wavelength to observe the molecular clouds. The CO 
transition (rotational quantum number 7=1-0) can be excited via collision with H2 
at 10 K. H2 is formed from Hi by reactions on dust surface and is then liberated into 
the interstellar space. There is no reliable tracer for a density range, approximately 
100-1000 cm^^. In this range, CO/H2 may not be formed yet and Hi may be opti- 
cally thick. It is therefore difficult to see the phase transition from Hi to H2, while 
future studies may be able to establish CI or CII transitions to trace such a density 
regime. 

The molecular clouds are clumpier than the Hi gas and the mass of a molecular 
cloud ranges from 100 M^, to 10^ M©. The large clouds having mass greater than 
10^ Mq are called giant molecular clouds (GMCs). The total number of GMCs in 
the Galaxy is estimated to be 3000. GMCs are nearly self-gravitating (e.g., Fukui 
and Kawamura 2010). They have long lifetime of a few 10 Myrs (Kawamura et al. 
2009) and are the formation site of stars of various mass. Ultraviolet radiation and 
stellar winds of high-mass stars formed in the GMCs ionize and eventually disperse 
GMCs. Such high-mass stars are the progenitors of SN explosions where the CR 
rays are accelerated. 
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3 Gamma-ray SNRs (supernova remnants) 

The interstellar medium (ISM) is crucial in producing high-energy radiation from 
supernova remnants. In particular, the protons in the ISM play a role as targets for 
CR protons to produce 7-rays. CR protons do not emit electromagnetic radiation, 
whereas a reaction between CR protons and ISM protons produces a neutral pion 
which decays into two 7-rays. This is the so called hadronic process of 7-ray produc- 
tion and provides a unique probe for the CR protons. Another reaction that produces 
7-rays is the leptonic process in which CR electrons energize low energy photons 
into 7-rays via the inverse Compton effect. Such CR electrons are observed directly 
via non-thermal synchrotron emission. One of the most important issues related to 
the 7-ray SNRs is if the 7-rays are produced by the hadronc process or the leptonic 
process. It is of vital importance to establish the hadronic-origin 7-rays in order to 
verify the origin of the CRs, a long standing puzzle since 19 1 2 when the CRs were 
discovered by V. Hess. 

If the hadronic process is working, we expect that the spatial distribution of the 
y-rays corresponds to that of the ISM protons. Some of the previous efforts how- 
ever did not find reasonable correspondence between the 7-rays and the ISM in the 
two most outstanding TeV 7-ray SNRs, RX J1713.7-3946 and RX J0852.0-4622 
(Aharonina et al. 2006b, 2007b). These previous works used only mm-wave CO 
data to estimate the ISM proton distribution. It may however be necessary to in- 
clude atomic protons also as the target protons if their density is high enough (Sec- 
tion 2). We employ the TeV 7-ray observations by H.E.S.S. at angular resolution 
of 0. 1 degrees, as well as the CO and Hi data obtained with the NANTEN2 4m 
telescope and ATCA (the Australia Telescope Compact Array) combined with the 
Parkes 64 m telescope, respectively. We present the results of comparisons between 
the TeV 7-rays and the ISM protons in the two SNRs, RX J1713.7-3946 and RX 
J0852. 0—4622, to demonstrate that the TeV 7-ray distribution well corresponds to 
the ISM proton distribution where both atomic and molecular protons are taken into 
account (Fukui et al. 2012, 2013). 



3.1 RX J 17 13.7 -3946 

RX J17 13.7 -3946 is the brightest TeV 7-ray SNR detected in the Galactic plane 
survey with H.E.S.S. (Aharonian et al. 2006a) and is a most promising candidate 
where the origin of the 7-rays may be established. The 7-rays higher than 10 TeV 
is detected here, suggesting that the CR protons responsible may have energy close 
to the knee, lO'^ eV, if the 7-rays are hadronic. The SNR was discovered in X-rays 
with ROSAT (Pfeffermann & Aschenbach 1996) and soon the X-rays are found to 
be non-thermal synchrotron emission with no thermal features (Koyama et al. 1997). 
TeV 7-rays were first detected by CANGAROO (Enomoto et al. 2002) and H.E.S.S. 
resolved the shell-like TeV 7-ray distribution (Aharonian et al. 2004, 2006b and 
2007a). Considerable work has been devoted to explore the 7-ray emission mecha- 
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nisms (Aharonian et al. 2006b; Porter et al. 2006; Katz & Waxman 2008; Berezhko 
& Volk 2008; Ellison & Vladimirov 2008; Tanaka et al. 2008; Morlino et al. 2009; 
Acero et al. 2009; Ellison et al. 2010; Patnaude et al. 2010; Zirakashvili & Aharo- 
nian 2010; Abdo et al. 2011; Fang et al. 2011). 

The molecular gas interacting with the SNR was discovered in the CO 7=1-0 
emission in Vlsr, the velocity with respect to the local standard of rest, around 
—7 km s^'. Figure 2 shows that the distribution of the CO emission in RX 
J1713.7-3946 well correlates with the X-ray distribution (Fukui et al. 2003); the 
northwestern rim of the X-ray coincides with the most prominent CO peaks (Fukui 
et al. 2003, Moriguchi et al. 2005, Fukui 2008, and Sano et al. 2010). This X-ray dis- 
tribution is now interpreted as caused by the interaction between the shock front and 
the molecular clouds as modeled in the magneto-hydrodynamical (MHD) numerical 
simulations by Inoue et al. (2012). The correlation provides a robust verification of 
the physical association of the CO clouds with the SNR shell. 

The distance of the SNR was determined to be 1 kpc by using the flat rotation 
curve of the Galaxy (Fukui et al. 2003). Studies of X-ray absorption suggested a 
similar distance 1 kpc under an assumption of uniform foreground gas distribution 
(Koyama et al. 1997), but the local bubble of Hi located by chance toward the SNR 
makes the X-ray absorption uncertain in estimating the distance (Slane et al.l999). 
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Fig. 2 Overlay map in Galactic coordinates showing a supernova remnant (SNR), RX 
J1713.7-3946, in gray scale [ROSAT PSPC X-ray Survey (Slane et al. 1999); from ROSAT archive 
database] and the intensity distribution of CO(7=1-0) emission in magenta contours. The inten- 
sity is derived by integrating the CO spectra from —1 1 to —3 km s^' , which is considered to be a 
velocity component interacting with the SNR. The lowest contour level and interval of CO are 4 
K km s-l (Fukui et al. 2003). 
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Fig. 3 (a) The H.E.S.S. TeV /-ray distribution of RX J1713.7— 3946 in smoothed excess counts 
above the cosmic-ray background (see Figure 2 of Aharonian et al. 2007a). Contours are plotted 
every 10 smoothed counts from 20 smoothed counts, (b) Averaged brightness temperature distri- 
bution of '^CO(7=1-0) emission in a velocity range of Vlsr = —20 km s^' to km s^' is shown 
in color (Fukui et al. 2003, Moriguchi et al. 2005). White contours show the H.E.S.S. TeV y-ray 
distribution and are plotted every 20 smoothed counts from 20 smoothed counts, (c) Averaged 
brightness temperature distribution of Hi emission obtained by ATCA and Parkes in a limited ve- 
locity range from Vlsr = —8 km s^ to —6 km s^' (McClure-Griffiths et al. 2005) is shown in 
color. White contours show the '^€0(7=1-0) brightness temperature integrated in the same veloc- 
ity range every 1.0 K km s^' (~3o') (Fukui et al. 2012). 



A subsequent careful analysis of the X-ray absorption also favors the smaller dis- 
tance (Cassam-Chenai et al. 2004). At 1 kpc the SNR has a radius of 9 pc and an 
age of 1600 yrs (Fukui et al. 2003; Wang et al. 1997) and the expanding shock 
front has a speed of 3000 km s^' (ZirakashviU & Aharonian 2010; Uchiyama et 
al. 2003, 2007). Most recently, Sano et al. (2010) showed that the SNR harbors the 
star forming dense clump peak C (Figure 3b) and argued that the X-ray intensity is 
significantly enhanced around the clump due to the shock interaction, reinforcing 
the association of the molecular gas. 

The molecular gas associated with the SNR opened a unique possibility to iden- 
tify target protons in the hadronic process. We expect the 7-ray distribution mimics 
that of the interstellar target protons if the hadronic process is working and if the 
CR distribution is uniform within the SNR. A detailed comparison between the ISM 
protons and the high-resolution 7-ray image of H.E.S.S. is therefore a useful test of 
the hadronic scenario. 

Figures 3a-c shows the distributions of TeV 7-rays, CO and Hi. The velocity 
range of the interacting gas with the SNR is later confirmed to be —20 km s^' 
- km s^' and is shown in Figure 3b (Fukui et al. 2012). We shall adopt this 
velocity range hereafter The TeV 7-rays show a well-defined SNR shell. The CO 
distribution shows a relatively good correspondence with the TeV 7-rays, whereas in 
some places the correspondence breaks. The most notable break is seen toward the 
southeastern rim (SB rim) of the 7-ray shell where no CO is seen. The Hi distribution 
shows that Hi is distributed over the entire SNR with a peak in the north, suggesting 
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Fig. 4 (a)The H.E.S.S. TeV y-ray distribution toward the SE cloud (Aharonian et al. 2007a). Red 
contours show averaged Hi brightness temperature distribution in a velocity range from —15 km 
s"^ to -5 km S-' (McClure-Griffiths et al. 2005). (b) The Hi and '^00(7=1-0) spectra at (/, b) 
=(347.55 degrees, —0.92 degrees). The shaded area shows an expected Hi profile. (Fukui et al. 
2012) 



that a significant amount of Hi gas is also associated with the SNR. The Hi intensity 
shows decrease toward the west and southeast in the SNR. The decrease in the west 
coincides with the CO and represents low-abundance cold Hi in the molecular cloud. 
The decrease in the southeast is different from that, because there is no CO in the 
SE rim. 

Figure 4a shows the Hi distribution overlayed on the TeV y-ray distribution to- 
ward the SE rim and Figure 4b shows a typical Hi profile and a CO profile in the 
SE rim at (Z, b) = (347.55 degrees, —0.92 degrees). The Hi contours show a fairly 
good correlation with the 7-rays, whereas the Hi intensity actually decreases toward 
the 7-ray enhancements. The Hi profile shows a broad dip in a velocity range from 
—20 to km s^'. We interpret that the dip is due to self-absorption by cool Hi as- 
sociated with the SNR. Such a broad dip is unusual for cold Hi dips in typical dark 
clouds which have generally a narrow Unewidth of a few km s^^. In the present 
case, we ascribe the broad dip as due to the large expanding motion of the Hi cav- 
ity wall. This Hi gas is not seen in CO, indicating that gas density is lower than 
10^ cm^^ and that the gas still contains fairly dense Hi in the order of 100 cm^^ 
as estimated below. The spin temperature T^ of cool Hi decreases with density due 
to the increased atomic line cooling and photon shielding as shown by numerical 
calculations (e.g.. Figure 2 in Goldsmith, Li & Krco 2007). This is just consistent 
with the behavior in Figure 4a; the decease of the Hi brightness indicates higher Hi 
density in the SE rim and the y-ray count increases with the density increase. By 
taking reasonable Ts ('-^40 K) as inferred from the minimum intensity of Hi at the 
dip bottom, we calculate the Hi column density responsible for the self-absorption. 
Here we adopt a conventional assumption that the background Hi emission is inter- 
polated by a straight Une (Figure 4b). We estimate that the Hi optical depth is close 
to 1 and calculate the absorbing Hi column density to be ^3 x 10^' cm^^ (see Figure 
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Fig. 5 (a) Distributions of column density of ISIVI protons Np estimated from '^CO(7=1-0), 
A'p(H2), (b) Hi emission with correction for tlie Hi self-absorption A'p(Hl) and (c) A'p(H2+Hl). 
All the datasets used here are smoothed to FWHIVI of TeV y-ray PSF. (d) TeV y-ray distribution. 
Contours are plotted every 50 smoothed counts from 20 smoothed counts (Fukui et al. 2012). 



5b), indicating that the SE rim has enhanced Hi density. As shown later in Figure 7, 
this interpretation of the Hi dip is consistent with the visual extinction which also 
shows enhanced ISM toward the SE rim. 

We present below some additional details on derivation of the ISM proton distri- 
bution. The molecular column density is calculated by using an X factor to convert 
the '^CO 7=1-0 intensity into H2 column density (equation 1) and the error is due 
to the 3(7 noise fluctuations in the CO observations. 



N{U2)=X-WC'iCO)) 



(1) 



where tiie X factor, N(ll2) I WC^CO){Yi km s"^), adopted is 2.0 x lO^" [cm^^ / 
(K km s~^)] (Bertsch et al. 1993). Then, the molecular proton density is given as A^p 

(H2) = 2A^(H2). 
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Fig. 6 (a) Distributions of column density of the total ISM protons WpCHi+Hl) in a velocity range 
from —20 km s^' to km s^' . Contours are the same as in Figure 1(a). (b) Azimuthal distributions 
of A'p(H2), A'p(Hl), A'p(H2+Hl) and TeV 7-ray smoothed counts per beam between the two elliptical 
rings shown in Figure 8(a). The proton column densities are averaged values between the rings 
(see text). Semi-major and semi-minor radii of the outer ring are 0.46 degrees and 0.42 degrees, 
respectively, and the radii of the inner ring are half of them. The same plots inside the inner ring 
are shown on the right side of Figure 8(b) (Fukui et al. 2012). 



The peak Hi brightness is usually around 80-100 K in the region except for the 
regions with Hi self-absorption. The brightness is consistent with warm Hi with 
typical spin temperature Ts ^125 K. The atomic proton column density is estimated 
by assuming that the 21 cm Hi line is optically thin as follows (Dickey & Lockman 
1990); 

A^p(Hl) = 1.823 X 10'*^ /rb^y (cm-2) (2) 

where T\, (K) and V (km s^') are the Hi brightness and velocity. The self- 
absorbing Hi column density is added to this in order to obtain the total atomic 
protons. The noise fluctuations in the Hi brightness are very small. The major un- 
certainty here comes from the uncertainty of T, of the cool Hi and the interpolation 
of the background Hi, but it is smaller than 10^' cm^^ and do not significantly affect 
the present estimates of A^p(H2H-Hi). 

The total ISM proton column density is finally given by the sum as follows; 



A^p(H2+Hl)=A?p(H2)+A^p(Hl) 



(3) 



The average H2 density of the CO cloud is around several 100 cm^^ and that of the 
SE rim is around 100 cm^^. 

Figure 5 summarizes the results. Figures 5a-d shows A^p(H2), A^p(Hl), A^p(H2H-Hl) 
and TeV 7-rays, respectively. First, we see that there is a significant amount of 
A^p(Hl) (column density is several times 10^' cm^^) over the SNR particularly to- 
ward the SE rim. The total ISM protons shows a significant difference by inclusion 
of the Hi and the ISM proton distribution shows a much better correspondence with 
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Fig. 7 (a)Ay distribution (Dobashi et al. 2005) is shown in color. TeV y-ray contours are the same 
as in Figure 1(a). (b) Distribution of column density of the total ISM protons estimated from both 
CO and Hi in a velocity range from Vlsr = ~20 to 10 km s^' . Here the Hi self-absorption is taken 
into account. TeV y-ray contours are the same as in Figure 1(a) (Fukui et al. 2012). 



the 7-ray distribution than the H2 alone. Figure 6 shows a comparison of the az- 
imuthal angular distribution between A'p and 7-rays and we see A^p(H2+Hl) shows a 
good correspondence with the TeV 7-rays. This is the first result which demonstrates 
a good spatial correspondence between TeV 7-rays and ISM protons at 0. 1 degree 
angular resolution and proves that a necessary condition for the hadronic scenario is 
fulfilled in the TeV 7-ray SNR. One of the natural next steps is to extend the method 
to the other SNRs as shown for RX J0852. 0—4622 in the next sub-section. 

An important implication of the Hi analysis is that the SE rim has enhanced col- 
umn density in Hi as in the TeV 7-rays, whereas the analysis of the Hi dip has some 
uncertainty in T^ and the Hi background interpolation. In order to make an inde- 
pendent check of the Hi column density, we can use other methods to measure the 
ISM protons. The visual extinction is one of such methods applicable at a distance 
of 1 kpc. Figure 7a shows that the distribution of the visual extinction is consistent 
with the Hi result in the SE rim, lending another support for the analysis of the Hi 
dip as self-absorption. Figure 7b shows for comparison the sum of the ISM proton 
distribution including both the foreground ISM protons and those in the SNR. More 
recently, a new analysis of the Suzaku X-ray data also shows a similar enhanced 
absorption column density toward the SE rim (Sano et al. 2013b). Considering all 
these, we are confident that the Ap(H2H-Hl) distribution in Figure 6a is reliable. 



3.2 RX J0852.0-4622 



RX J0852.0-4622 is another prominent young TeV 7-ray SNR imaged by H.E.S.S. 
Aschenbach et al. (1998) discovered RX J0852.0-4622 as a hard X-ray SNR 
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Fig. 8 (a) The H.E.S.S. TeV y-ray distribution of RX J0852.0-4622 in smoothed excess counts 
(see Figure 1 of Aharonian et al. 2007b). Contours are plotted every 10 smoothed counts from 50 
smoothed counts, (b) Integrated intensity distribution of '^CO(/=1-0) emission in a velocity range 
of ^LSR = 24 km s^' to 30 km s^' is shown in color (Moriguchi et al. 2001). White contours show 
the H.E.S.S. TeV y-ray distribution and are plotted the same as Figure 8 (a). 



seen in projection against the Vela SNR in ROSAT All-Sky Survey image. RX 
J0852. 0—4622 and RX J1713.7— 3946 share similar properties; they are both young 
with ages of 1600-2000 yrs and show non-thermal synchrotron shell-like X-ray 
emission with no thermal features in addition to the strong shell-like TeV 7-rays 
(Figure 8a). The apparent large diameter of RX J0852. 0—4622, about 2 degrees, is 
most favorable to test a spatial correspondence between the 7-rays and the ISM at 
a 0.12 degrees angular resolution (FWHM) of H.E.S.S.. There are two more simi- 
lar shell-like TeV 7-ray SNRs, RCW 86 and HESS J1731-347, but they are small 
with a size of about 0.5 degrees (Aharonian et al. 2009; H.E.S.S. Collaboration et al. 
201 1), making it difficult to test the spatial correspondence. 

The age and distance of RX J0852. 0—4622 have been the subject of debate in 
the hterature. Slane et al. (2001a, b) argued that RX J0852. 0-4622 is physically 
associated with the GMC, the Vela Molecular Ridge (VMR, May et al. 1988; Ya- 
maguchi et al. 1999). The distance of the VMR is estimated to be 700 ± 200 pc 
(Liseau et al. 1992), whereas it is not yet established if the VMR is physically con- 
nected to RX J0852.0-4622 (see e.g., Pannuti et al. 2010). Recently, Katsuda et 
al. (2008) estimated an expansion rate of RX J0852. 0—4622 based on two observa- 
tions separated by 6.5 years made with XMM-Newton toward the northwestern rim 
of RX J0852.0-4622, and derived an age of 1.7-4.3 x 10^ years for the SNR and a 
distance of 750 pc. We will adopt this distance to RX J0852.0-4622. 

In the region of RX J0852. 0—4622, the CO clouds are significantly less than in 
the region of RX J1713.7-3946. This is not surprising since RX J0852.0-4622 is 
far from the Galactic center where the ISM density is generally lower than in the 
region of RX J1713.7— 3946. Among all the CO features in the direction, we find 
that the CO clouds at 25 km s^' show clear association with the southwestern rim 
of the SNR (Figure 8b). Based on this correspondence we have searched for asso- 
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Fig. 9 (a) The column density of the ISM proton distribution A'p(H2+Hl) and the contours of TeV 
/-rays. The outer and inner circles are centered at (/, b) = (266.28 degrees, —1.24 degrees) with 
radii of 0.92 degrees and 0.63 degrees, respectively, (b) The azimuthal angles are indicated. The 
azimuthal distribution of the TeV /-rays and the ISM protons. The /-ray counts and the column 
density of the ISM protons in the two circles in Figure 9 (a) are shown at every 30 degrees. The 
black dots show the TeV /-ray counts and the blue box the molecular protons A''p(H2), the green 
box the atomic protons A'p(Hl), and the red triangle the total ISM protons Ai'pfH^-FHl). The right- 
most point shows the average within the innercircle (Fukui et al. 2013). 



dated CO and Hi gas with the SNR shell. Figure 9a shows the total ISM proton 
distribution in RX J0852. 0—4622 derived by combining the associated CO and Hi 
distributions. The velocity range is from to 50 km s^' centered at 25 km s^', al- 
though the Galactic rotation suggests lower velocity less than 10 km s^'. We ascribe 
this velocity deviation as due to the expansion of a few Hi supersheUs; the SNR RX 
J0852. 0—4622 is located toward an overlapping region among these Hi shells and 
the associated ISM has significant deviation by 50 km s^' at maximum from the 
pure Galactic rotation. The total associated Hi mass is 10^ Mq while that of molec- 
ular gas is 10^ M,7,. So, the associated ISM is dominated by the atomic gas in RX 
J0852.0-4622. On the other hand, the associated ISM mass in RX J1713.7-3946 is 
10"* Mq both in H2 and Hi, respectively. RX J1713. 7-3946 is therefore associated 
with 10 times more molecular gas than RX J0852. 0—4622. 

We made a Gaussian fitting to the radial TeV 7-ray distribution and determined 
the best-fit parameters Uke the radius and width of the shell (Figure 9a). Figure 9b 
shows the azimuthal TeV 7-ray distribution (Fukui et al. 2013). This is obtained 
by averaging the 7-ray count between the two rings at the 1/3 level of the peak of 
the Gaussian shell, 0.61 degrees and 0.92 degrees in the radial distribution from 
the SNR center The two rings include the major part of the 7-ray shell and Fig- 
ure 9b represents well the 7-ray azimuthal distribution. We find the correspondence 
between the total ISM and the TeV 7-rays is reasonably good. This is the second 
case that shows a good correspondence between the ISM protons and TeV 7-rays. 
Toward the angle 15 degrees in Figure 9b there is some excess in the 7-rays. This 
is due to a PWN located by chance toward this part of the shell, which and is not 
associated with the SNR. 
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To summarize this section, we have shown that the ISM proton distribution 
A^p(H2+Hl) shows a good spatial correspondence with the TeV 7-ray distribution 
in the two outstanding TeV 7-ray SNRs. This is a step forward in establishing the 
hadronic origin of 7-rays, whereas the results alone do not reject the leptonic sce- 
nario. What remains to be explored is 1) if a leptonic model for the 7-ray production 
can also explain the observed TeV 7-ray distribution, and 2) if such a hadronic model 
is consistent with the non-thermal X-ray distribution by the CR electrons. 



4 Interaction between the ISM and the SNR shock, and the y-ray 
production 

The comparisons between 7-rays and the ISM in the two SNRs in Section 3 have 
led to a picture that the shell-like clumpy ISM consisting of atomic and molecu- 
lar gas acts as the target for CR protons to produce the hadronic 7-rays. We shall 
here discuss into more detail the interaction between the SN shock front and the 
ISM in RX J1713.7— 3946. We presume that a similar argument is applicable to RX 
J0852. 0—4622 by considering its properties similar to RX J1713.7— 3946 (Section 

3). 

Figure 10a shows a comparison between the non-thermal X-rays and the CO 
distribution in the northwestern (NW) rim of RX J1713.7-3946 (Sano et al. 2010). 
In the X-ray distribution, the NW rim is particularly intense and this is where the 
densest CO clumps are located. Figure 10b shows a close-up view toward peak C 
where a protostar with bipolar outflow is located in the center Figure 10 shows that 
the X-rays are correlated with CO at pc scale but is anti-correlated at sub-pc scale 
as shown by that peak C is located toward the local minimum of the X-rays. This 
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Fig. 10 Suzaku XIS (XIS 0+2+3) mosaic image of RX J1713.7-3946 in the energy bands 1-5 keV 
in color scale and (a) '2cO(7=2-l) (b) ^^CO{J=4-3) intensity contours obtained by NANTEN2 
(Sano et al. 2010 Figure 5). The area enclosed by the blackbox in (a) is shown enlarged in Figure 
9 (b). Open crosses show the positions of the IRAS point sources. 
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behavior is seen in every CO peak within this SNR (Sano et al. 2013a), offering a 
hint on the details of the interaction. 

The initial ISM distribution prior to the SN explosion consists of a wind- 
evacuated low-density cavity of ^9 pc radius and the cavity wall consists of around 
ten dense molecular clumps, the cool and dense Hi gas as in the SE rim, and the 
warm Hi (Figures 10 and 11). The interaction with the SN shock waves is taking 
place in the last several 100 yrs as estimated from a ratio of ~3 pc / 3000 km s^'. 
In such a short time scale it is impossible that the molecular clumps are pushed ap- 
preciably in space and it is likely that the neutral gas distribution, in particular the 
molecular distribution, is not significantly affected by the shock interaction. Some 
of the dense CO clumps located in the inner part of the shell including peak C are 
likely overtaken by the shock (Sano et al. 2010). In addition, the stellar winds likely 
have stripped off the lower density envelope of the CO clumps, making a steep 
density gradient at the clump surface. Such a gradient is in fact observed toward 
peak C by using the sub-mm transitions of CO to trace the radial molecular density 
distribution (Sano et al. 2010). 

The SN shock front having 3000 km s^' expansion velocity interacts with the 
ambient highly clumpy ISM distribution. Many of the previous works on the shock 
interaction in the SNR assume that the ISM is uniform in density for simplicity 
(e.g., Ellison et al. 2010). If we assume such uniform ISM, in order to produce the 
observed 7-rays the ISM density has to be increased to at least 10 cm^^. Then, the 
uniform model predicts that strong thermal X-rays will be emitted by the heating of 
the ISM in the shock waves. This creates a serious discrepancy with observations 
since no thermal X-rays are detected. This was one of the major arguments against 
the hadronic scenario in the RX J1713.7-3946 (e.g., Ellison et al. 2010). While 
the uniform ISM may be possible as a first approximation, more accurate properties 
of the interaction strongly depend on the ISM density inhomogeneity. We should 
therefore adopt more realistic highly inhomogeneous ISM distribution as observed 
in CO (Figure 3). 

Inoue et al. (2012) carried out MHD numerical simulations by considering the 
clumpy ISM distribution consistent with the observations in RX J1713.7— 3946 and 
discussed the basic aspects of the interaction and the y-ray and X-ray properties. The 
interaction between the shock and the CO clumps creates turbulence as simulations 
indicate (Figure 12). This turbulence amplifies the magnetic field around the clumps 
and enhances the synchrotron X-rays. The magnetic field is in the order of 10-100 
JJ.G on the average but may be enhanced to the mG order (Uchiyama et al. 2007). 
The shock waves cannot penetrate into the dense molecular clumps and cannot heat 
them up. This explains the lack of thermal X-rays which happens usually in the 
other SNRs of lower-density environments. The CR electrons cannot penetrate into 
the dense clumps since its penetration depth is only in the order of ~0. 1 pc in a 
given timescale ~1000 yrs (Fukui et al. 2012). As a result, the X-rays are enhanced 
around the CO clump but are depressed inside the clumps. This offers an explanation 
of the pc scale correlation and sub-pc scale anti-correlation between CO and X- 
rays (Figure 10). The inside of the cavity is as a whole of very low density due 
to the stellar winds evacuation, where particle acceleration is efficiently made via 
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diffusive shock acceleration (DSA) to the energy close to the knee (e.g.. Bell 1978; 
Blandford & Ostriker 1978 ). CR protons have larger penetration lengths of ^1 pc 
than electrons and can interact with the ISM protons inside the molecular clumps. 
Such proton-proton reaction explains the y-ray production and is consistent with the 
spatial correspondence between the 7-rays and the ISM protons shown in Figures 
6 and 9. The inhomogeneous ISM therefore solves the difficulties in the uniform 
ISM models by allowing the y-ray production in the dense clumps as well as the 
suppression of the thermal X-rays. 

In order to have a deeper insight into the 7-ray production, we need to understand 
the 7-ray energy spectrum. It is usually assumed that CR protons can interact with 
all the target ISM protons in a SNR, but, to be realistic, penetration of CR protons 
into the molecular gas is energy dependent for a high density regime like ■~10'* cm^^ 
or more. The high-energy CR protons can penetrate into the dense clumps, whereas 
the low-energy protons cannot (Gabici et al. 2007, Zirakashvili & Aharonian 2010). 
This energy dependent penetration makes the lower energy protons interact only 
with the lower density ISM protons. This reflects in the 7-ray energy spectrum and 
a harder 7-ray spectrum is expected than in the uniform ISM case. 
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Fig. 11 Schematic picture of the shock-cloud interaction model (Inoue et al. 2012 Figure 10). The 
primary forward shock wave propagates through the cloudy wind bubble, where particle accelera- 
tion operates by the diffusive shock acceleration (DSA). The shock waves in the clouds are stalled, 
which suppresses thermal X-ray line emission. Shock-cloud interactions induce turbulent eddies, 
which amplifies the magnetic field (~1 mG) caused by the turbulent dynamo effect and the reverse 
shocks. Then the synchrotron emissions enhanced ai'ound the clouds where magnetic field strength 
is ~ 1 mG around shocked clouds. Also, hadronic y-rays are emitted from dense clouds illuminated 
by accelerated protons whose photon index can be p— ill = 1.5 for p = 2. 
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(a) time = 750 yr 




Fig. 12 Results of the numerical simulations between the SN Shockwaves on the ISM (Inoue et al. 
2012). The perpendicular shock case at f = 750 yr after the shock injection (Inoue et al. 2012). (a): 
number density volume rendering. Regions in green and blue indicate the density n ~10 cm^-' and 
n > 30 cm^^^, respectively, and the regions in warm colors show the shocked diffuse gas with n < 
4 cm^ . Magnetic field lines are represented as gray lines, (b): two-dimensional number density 
slice at z = 0.0 pc. (c): close up of the solid box in (b). (d): volume rendering of magnetic field 
strength, (e): slice of magnetic field strength at z = 0.0 pc. (f): close up of the solid box in (e). 



Figures 13 shows the 7-ray spectrum of RX J1713.7— 3946 covering the energy 
range 10^-10** MeV (Abdo et al. 2011). This is the combined data of H.E.S.S. and 
Fermi satelhte (Aharonian et al. 2007a; Abdo et al. 201 1). RX J1713.7-3946 shows 
a hard spectrum which is not compatible with the conventional hadronic y-rays and 
these authors claimed that the spectrum supports the leptonic scenario. The density- 
and energy-dependence of the hadronic interaction however indicates that such a 
hard spectrum is also explained well in the hadronic scenario with highly inhomoge- 
neous ISM distribution as argued into detail by Inoue et al. (2012). Figure 14 shows 
the 7-ray spectrum in the same energy band with Figure 13 for RX J0852. 0—4622 
(Tanaka et al. 2011). This shows a softer spectrum than in Figure 13 and may be 
better explained by the hadronic scenario than the leptonic scenario. In connection 
with the ISM distribution, we note that RX J0852. 0-4622 has much less CO clouds 
than in RX J1713.7— 3946. CR protons will therefore interact mainly with the Hi 
gas whose mass is 10 times more than the molecular mass, offering an explanation 
of the softer spectrum by the hadronic scenario. These two cases suggest that the 



Molecular and atomic gas in the young TeV y-ray SNRs 



17 



ISM distribution is crucial in interpreting the 7-ray properties and, in particular, in 
exploring the hadronic 7-ray s. 

The identification of target ISM protons in the hadronic scenario allows us to esti- 
mate the total CR proton energy Wp in these SNRs. By using the relationships which 
connect the 7-ray counts and the ISM density (Aharonian et al. 2006b, 2007a), we 
estimate that Wp is 10"*^ ergs in the two SNRs. This energy corresponds to 0.1% of 
the total kinetic energy of a SNR and seems to be 10 times less than what is required 
to explain the CR energy budget of the whole Galaxy. There are a few possibilities 
to reconcile this discrepancy. The CR energy may yet increase in time with the evo- 
lution of a SNR and the total CR energy may become more than estimated above. 
In W44, a middle-aged SNR, the lower limit of Wp is estimated to be 10"*^ erg, 10 
times larger than the present value (Giuliani et al. 2011, Yoshiike et al. 2012). It is 
also to be considered that the ISM volume filling factor in a SNR shell may not be 
100%. Particularly, in RX J0852.0-4622 the eastern half of the SNR does not have 
much ISM, whereas we expect CR protons have a similar number density over the 
SNR. This suggests that some portion of the CR protons may not be detected in the 
7-rays, leading to an under-estimate of Wp by some factor. 

In order to better understand the origin of 7-rays in SNRs, we need to have 
a larger sample of SNRs. A direction is to extend an analysis of the ISM to the 
middle-aged SNRs like W44, and the other to study smaller TeV 7-ray SNRs like 
HESS J1731-347 (H.E.S.S.Collaboration et al. 2011). The latter should be possi- 
ble with the future instrument CTA at a higher angular resolution of 1 arcmin with 
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Fig. 13 Energy spectrum and numerical models of RX J 17 13. 7— 3946 in GeV to TeV y-ray spectra 
taken by the Fermi Large Area Telescope and H.E.S.S. (Abdo et al. 2011 Figure. 3). Some solid 
lines indicate the numerical model of the hadronic model (;r"-decay spectrum: left) and the lep- 
tonic model (inverse Compton scattering spectrum: right). The green region shows the uncertainty 
band obtained from maximum likelihood fit of the spectrum assuming a power law between 500 
MeV and 400 GeV for the default model of the region. The gray region depicts the systematic 
uncertainty of this fit obtained by variation of the background and source models. The black error 
bars con'cspond to independent fits of the flux of RX J1713.7— 3946 in the respective energy bands. 
Upper limits are set at a 95% confidence level. 
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Fig. 14 Energy spectrum and numerical models of RX J0852.0— 4622 in GeV to TeV y-ray spectra 
taken by the Fermi Large Area Telescope (red) and H.E.S.S. (black) (Tanaka et al. 2011 Figure.2 
and 3). For the Fermi-LAT points in left panel, the vertical red lines and the black caps indi- 
cate statistical and systematic errors, respectively. The dotted line shows the best-fit power-law 
obtained from the maximum likelihood fit for the entire 1-300 GeV band (The butterfly shape 
is the 68% confidence region). The dashed curve is the ;r''-decay spectrum (hadronic model) by 
Berezhko et al. (2009). Broadband SED of RX J0852.0-4622 with (a) the hadronic model and (b) 
the leptonic model. The radio data points and blue line indicate the integrated fluxes of the SNR 
determined based on the 64 m Parkes radio telescope (Duncan & Green 2000) and the the X-ray 
flux using the ASCA data-set (Aharonian et al. 2007b). The solid, dashed, and dot-dashed lines 
represent contributions from ;r''-decays, inverse Corapton scattering, and synchrotron radiation, 
respectively. 



higher sensitivity. This will definitely help us to increase the number of samples 
significantly and provide a better understanding of the 7-ray origin. 



5 Dark gas 



An important topic related to the 7-rays originating from the ISM protons is dark 
gas. The term dark gas is used to represent unseen gas either in Hi or CO but is 
detectable as y-ray excess (Grenier et al. 2005). The first observational study of dark 
gas was made by using the EGRET data toward the local ISM and showed that there 
is y-ray excess, at galactic latitude 5 degrees < \b\ < 80 degrees, beyond the level 
which the gas detected by the CO and Hi transitions is able to explain. Recently, 
the Planck satellite provided sensitive mm/sub-mm distributions of the same region 
and showed that 15% of the mm/sub-mm dust emission is not explained either by 
the CO or Hi gas (Planck Collaboration et al. 201 lb). This is also called dark gas. 
In addition, Planck satellite presents similar dark gas components in the Magellanic 
Clouds (Planck Collaboration et al. 2011a). It is often suggested that dark gas may 
corresponds to molecular gas with no detectable CO (Planck Collaboration et al. 
2011b), while another possibility discussed is opaque Hi gas whose intensity is 
saturated by moderately large optical depth. 
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Fig. 15 ATi, and T variations depending on A'p(Hl) shown by Torii et al. (2012). Here AV of 10 
km s^' is assumed. The filled area shows a Np range of the dark gas (Grenier et al. 2005; Planck 
Collaboration et al. 2011b). 



The results of the 7-ray SNR in Section 3 suggest that cool Hi gas having 
large optical depth around 1 is responsible for producing part of the 7-rays via the 
hadronic process. We here discuss that dark gas may correspond to such cool and 
dense Hi gas with relatively large optical depth. Figure 15 shows the Hi intensity 
as a function of Hi column density for two spin temperatures 7^= 100 K and 40 K. 
The assumption here is that the Hi linewidth is 10 km s^'. The plot shows that the 
Hi intensity becomes saturated in the Hi column density range from 10^' cm^^ to 
10^^ cm^^. For T^ = 40 K, the optical depth of Hi is greater than 1 and the typical 
Hi density is around 100 cm^^. These Hi column densities correspond to that of the 
dark gas observed by Planck and we suggest that the purely atomic dense gas is a 
viable alternative to CO-free molecular gas to explain the dark gas. 

We discuss on the time scale of H2 formation in order to shed another light on 
the dark gas composition. The H2 molecules are formed on dust surface with the 
formation time scale ffonn given by the equation (4) (HoUenbach & Salpeter 1971). 



fform '^ 10 



100 



n (cm ^) 



(yr) 



(4) 



For density of n= 100 cm^ , the H2 formation timescale is 10 Myrs. This is similar 
to the typical evolutionary timescale of GMCs. On the other hand, the local gas 
at high galactic latitude where the dark gas is identified is all low column density 
clouds whose Hi column density is in the order of 10^' cm^^. They are surrounding 
well known local molecular clouds like Ophiuchus and Lupus. The crossing time 
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scale of such ensembles of local clouds, most likely not gravitationally bound, is 
estimated to be ^1 Myrs as derived from a ratio 10 pc / 10 km s^' (e.g., see Figure 
4 in Tachihara et al., 2001). So, the time scale of the surrounding gas is too short 
to convert Hi into H2. This argues for that the dark gas is dominated by dense Hi 
but not by H2 without CO. On the other hand, it is conceivable that the envelops of 
the GMCs in the Galactic plane may have lifetime long enough to form H2, and H2 
envelope without CO may be a viable alternative. Such GMCs are not well sampled 
in the Galaxy due to heavy contamination in the Galactic plane but the GMCs in the 
LMC (Kawamura et al. 2009) may provide an observable case where the molecular 
dark gas without CO is dominant. Future careful studies of the dust components will 
help to clarify the real contents of the dark gas. 



6 Summary 

From a view point of high-energy astrophysics, in particular with an emphasis of 
the origin of 7-rays, we present properties of neutral gas in the interstellar medium 
and its connection with the y-ray SNRs. 

We first review the gas properties of both atomic and molecular forms as probed 
by the 2.6 mm CO and 21 cm Hi transitions. CO traces dense molecular gas and Hi 
traces lower density atomic gas, whereas either CO or Hi may not be a good tracer 
in an intermediate density regime roughly between 100 cm^^^ and 1000 cm^^. 

We then present analyses of the two most remarkable TeV 7-ray SNRs RX 
J1713.7-3946 and RX J0852. 0-4622 and demonstrate that the ISM proton dis- 
tribution shows a good spatial correspondence with the TeV 7-ray distribution. We 
argue that this correspondence provides a support for the hadronic scenario of 7-ray 
production. What remains to be explored is if a leptonic model consistent with the 
non-thermal X-ray distribution can also explain the observed TeV 7-ray distribution. 
We describe theoretical numerical simulations of the interaction of the SNR shock 
waves with the clumpy ISM and discuss that such realistic simulations can explain 
the observations reasonably well, whereas the uniform ISM models are not viable. 

Then, we present a discussion on the dark gas observed in the 7-rays and mm/sub- 
mm dust emission. We argue that, for the local component of the dark gas in the 
Galaxy, dense Hi gas as observed in RX J1713.7— 3946 is a good candidate for 
the dark gas. The local gas is dominated by low mass clouds and their time scale 
may be too short to convert Hi into H2. On the other hand, the GMCs in the LMC 
have larger time scale like 10 Myrs, long enough for H2 formation, and therefore 
molecular dominant dark gas is a viable model on the GMC scale. 

Finally, we conclude that the ISM plays an essential role in producing the 7-rays 
and X-rays in the SNR, and the future instruments like CTA, combined with the 
ISM studies, will achieve tremendous progress in the 7-ray astrophysics of the SNR 
and the other high-energy objects. 
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